; . 
• 1 J, 


AFWL-TR-72-19 


SOME FUNDAMENTAL ASPECTS OF 
NUCLEAR WEAPONS 


Henry F. Cooper, Jr. 


Harold L. Brode 


Gerald G. Leigh 


TECHNICAL REPORT NO. AFWL-TR-72-19 


March 1972 
r---. 
~--------_____ 
_ 
J---.-- - 
~ 
"" 


I I . 
~ 
__ 
~========r-'- 
AIR FORCE WEAPONS LABORATORY 


Air Force Systems Command 


Kirtland Air Force Base 


New Mexico 


AFWL-TR- 
72-19 


DOC 


Approved for public release; distribution unlimited. 


;'; 


AIR FORCE WEAPONS LABORATORY 


Air Force Systems Command 
Kirtland Air Force Base 
New Mexico 
87117 


When US Government drawings, specifications, or o.ther data are used for 
any purpose other than a definitely related Government procurement operation, 
the Government thereby incurs no responsibility nor any obligation whatsoever, 
and the fact that the Government may have formulated, furnished, or in any way 
supplied the said drawings, specifications, or other data. is not to be 
regarded by implication or otherwise, as in any.manner licensing the holder 
or any other person or corporation, or conveying any rights or permission to 
manufacture, use, or sell any patented invention that may in any way be 
related thereto. 


DO NOT RETURN THIS COPY. 
RETAIN OR DESTROY. 


'. 


UNCLASSIFIED 


Se 
CI 
If 
cunty 
asS Lcatlon 
OOCUMEIH CONTROL OAT A " R & 0 


(Security c-rassi!icatron 01 "Ue, body ot abstract and inde"ins annotation must be entered when the overall report Is clesslf/ed) 


I. ORIGINATING ACTIVITY (Corporate author) 
2a. REPORT SEC:URITy CL.ASSIFICATION 


Air Force Weapons LRboratory (DEV) 
UNCLASSIFIED 


Kirtland Air Force Base, New Mexico 
87117 
2b. GROUP 


3, REPORT TITLE 


SOME FUNDAMENTAL ASPECTS OF NUCLEAR WEAPONS 


4, DESCRIPTIVE NOTES (Type of report and inclusive dates) 


June 1971-AulZust 1971 


~. AU THORIS) (Firs! name, middle Initial, last name) 


Henry F. Cooper, Jr. ; Gerald G. Leigh, Major, USAF, Air Force Weapons Laboratory; 
Harold L. Brode, Research and Development Associates, Santa Monica, California 


e· REPORT DATE 
7a. TOTAL NO. OF PAGES 
rb. NO. OF REFS 


March 1972 
46 
18 
Ba. CONTRACT OR GRANT NO. 
9a. ORIGINATOR"S REPORT NUMBERIS) 


b. PROJECT NO. 
5710 
AFWL-TR-72-l9 


c. 
ab. OTHER REPORT NO<S) (Any other numbers thal'may be assIgned 
this reporl) 


d. 


10. DISTRIBUTION STATEMENT 


Approved for public release; distribution unlimited. 


I" SUPPLEMENTARY NOTES 
12. SPONSORING MILITARy"A'CTIVITY 


AFWL (DEV) 
r 


Kirtland AFB, NM 
87117 


~"'BSTRACT 
(Distribution Limitation Statement A) 


; Nuclear weapons effects associated with a surface burst geometry, particularly on 
a rock medium, are summarized. 
Particular emphasis is placed on the close-in blast 
and shock phenomenology that usually is the bounding factor limiting the potential 
hardness of protective structures. 
Consideration is given to the coupling of energy 
into the air and earth in the vicinity of ground zero, the close-in airblast and 
ground shock, and the subsequent environment usually occurring at ranges of 
strategic interest. ( ) IF-- 


. -------. 
Details of illustrations In . 
this document may be. better 
studied on mICrofiche 


UNCLASSIFIED 


Security ClassHication 


, ,. 
LINK A 
LIN K B 
LINK C 
KEY WORDS 
ROLE 
WT 
ROLE 
WT 
ROLE 
WT 


Nuclear weapons effects 
Blast and shock airblast 
Ground shock 
Ground motions 
Finite difference methods 


.- 


, .e"' 


, 


UNCLASSIFIED 


Security Classification 


AFWL-TR-72-19 


SOME FUNDAMENTAL ASPECTS OF NUCLEAR WEAPONS 


Henry F. Cooper, Jr. 


"Harold L. Brode 
Gerald'G. Leigh 


TECHNICAL REPORT NO. AFWL-TR-72-l9 


Approved for public release; distribution unli~ited. 


AFWL-TR-72-l9 


FOREWORD 


This research was performed under Program Element 6ll02H, -Projec~,57l0. 


,';,' ' 


Inclusive dates of research were June 1971 through August 1971. 
The report 
was submitted 18 January 1972 by the Air Force Weapons Laboratory Project 
Officer, Dr. Henry F. Cooper, Jr. (DEV). 


This report is an edited version of a paper presented at a NATO Defense 
Group Seminar in Bodo, Norway, 15 September 1971, by Dr. Henry F. Cooper, Jr., 
and Major Gerald G. Leigh of the Air Force Weapons Laboratory, Kirtland AFB, 
New Mexico, and Dr. Harold L. Brode of Research and Development Associates, 
Santa Monica, California. 


This technical report has been reviewed and is approved . 


.4~Lt/ 
HENRY F. COOPER, JR. 
Project Officer 


vI!,A;J t~ 


GERALD G. LEIGH 
Maj or,' USAF 
Chief, Facilities Survivability 
Branch 


:ff~?5~ 
WILLIAM B. LIDDICOET 
Colonel, USAF 
Chief, Civil Engineering Research 
Division 


ii 


AFWL-TR-72-l9 


ABSTRACT 


(Distribution Limitation Statement A) 


Nuclear weapons effects associated with a surface-burst geometry, particularly 
on a rock medium, are summarized. 
Particu~ar emphasis is placed on the close-in 
blast and shock phenomenology that usually is the bounding factor limiting the 
potential hardness of protective structures. 
Consideration is given to the 
. 


coupling of energy into the air and earth in the vicinity of ground zero, the 
close-in airblast and ground shock, and the subsequent environment usually 
occurring at ranges of strategic interest. 


iii/iv 


/ 


AFWL-TR-72-l9 


Section 


I 


II 


III 


IV 


V, 


CONTENTS 


INTRODUCTION 


THEORETICAL PREDICTION TECHNIQUES 


AIRBLAST AND ASSOCIATED PHENOMENA 


Airblast-Induced Ground Motions 


Debris and Fallout 


FREE-FIELD GROUND MOTIONS FROM SURFACE AND 
NEAR-SURFACE BURSTS 


Initial Energy Coupling 


Direct-Induced Ground Motions 


Crater-Induced Ground Motions 


SUMMARY AND CONCLUDING COMMENTS 


REFERENCES 


Preceding Page Blank 
v 


1 


3 


7 


9 


14 


18 


18 


21 


33 


36 


38 


/ 


AFWL-TR-72-l9 


Figure 


1 


2 


3 


4 


ILLUSTRATIONS 


Comparison of Peak Stress (0.5 kilobar) for a 
Megaton Burst on Rock 


Shock Radius (Rs) and Arrival Time (ts ) versus 
Peak Overpressure for a l-MT Surface Burst 


Temperature versus Time at High Peak Overpressures 
(l-MT surface burst) 


Shock-\-lave Peak Pressures (l-MT surface burst) 


I 
5 
Air-Shock-Induced Ground Motion Wave Fronts for 


Pea~ Overpressur¢s of 10,000, 1,000, 300, and 100 
psi and for Seis.fuic Velocities of 2,500 ft/secand 
5,000 ft/sec 


6 
Late Fireball Turbulence 


7· 
Initial Conditions for Surface-Burst Cratering 
Motion 


8 
Peak-Pressure Isobars for Surface 790-KT Sour~e 


9 
Peak-Pressure Isobars for I-Meter Burial of a 
790-KT Source 


, 
10 
Peak-Pressure Isobars ifor 8-Meter Burial of a 
790-KT Source 


11 
Comparison of Peak-Particle Velocity from Contained 
Bursts in Hard and Soft Rock 


12 
Phenomenological Regions for the Surface-Burst 
Geometry 


13 
Comparison of Peak-Particle Displacement from 
Contained Bursts in Hard and Soft Rock 


14 
Mineral Rock Particle Velocity at 70-Foot Range and 
2-Foot -Depth 


15 
Free-Surface Permanent Displacements (Mineral Rock) 


16 
Crater-Induced Surface Motion on Distant Plain 6 


17 
Estimates of Grater-Induced Displacements from 
Contact Bursts on Hard Rock 


vi 


6 


8 


10 


10 


12 


·16 


20 


22 


" 


23 


24 


25 


27 


30 


31 


32 


34 


35 


AFWL-TR-72-l9 


. SECTION I 


INTRODUCTION 


Most nuclear weapons effects that are pertinent to evaluating the hardness 


of various strategic systems have been summarized in the official unclassified 


source, The Effects of Nuclear Weapon8~ edited by Glasstone (Ref. 1). 
More 


recently, while emphasizing close-in nuclear explosion physice, Brode (Ref. 2) 


.-; 


provided a series of formulae for estimating various transient phenomena. 
These 


two complementary technical papers provide an excellent review of most nuclear 


weapons effects pertinent to the field of protective construction. 


This report will not dwell on the phenomenology contained in the above- 


mentioned references. 
Rather, it concentrates on more recently obtained infor- 


mation not contained in references 1 and 2. 
In particular, it reviews our cur- 


rent understanding of ground shock phenomena in soil and rock media. 
Other 


·effects are cursorily summarized in the context of near-surface explosions. 


The state-of-the-art understanding of ground motions 10 years ago--when the 


criteria were developed for many currently existing strategic structures--was 


such that the basic ground motion phenomena were known in only very general 


terms (Ref. 3). 
In addition, the dynamic interaction of soil with a structure 


had not been thoroughly studied and only rough approximations were provided to 


define the stnlctural landing and motion. 


Although the current understanding of these phenomena is still·incomplet~. 


we have progressed markedly in the last few years. 
The new understanding has 


come from theoretical and experimental programs aimed at increasing our knowledge 


of the basic phenomena and providing the tools with which to simulate such 


effects. 
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The scope of this report is limited to some nuclear weapons effects relevant 


to the hardness of strategic systems located on or near the earth's surface. 


We concentrate on surface and near-surface burst geometries, particularly for a 


rock medium. 
The effects considered include initial coupling phenomena, crater- 


ing, airblast, and ground motions. ' Initial coupling and airblast phenomena are 


fairly independent of geologic conditions,· but cratering and ground motions are 


strongly dependent on in-situ soil and rock properties. 
Thus, this report also 


emphasizes te~t results that clarify our understanding of the in-situ response 


of geologic media, particularly rock. 
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SECTION II 


THEORETICAL PREDICTION TECHNIQUES 


Before beginning our description of t~e basic phenomena, we will first 


review some aspects of the theoretical tools that play a major role in our 


current research programs. 


In recent years finite difference techniques pioneered by personnel at the 


Los Alamos Scientific 'Laboratories (LASL) and the Lawrence Livermore, Laboratories 


(LLL) have been further developed and applied to study airb1ast, cratering, ground 


motion, and other nuclear weapons effects. 
These theoretical procedures are 


based on first principles and treat radiative phenomena (transport and diffusion) 


coupled with material continuum response. 
The first major calculation of the 


surface-burst geometry (Ref. 4) considered the earth's respons~ under high 


pressure and treated the soil: as a hydrodynamic, fluid. 
In the past several 


years the theoret,ica1 techniques have been extended to consider material strength 


as well as hydrodynamic behavior. 


Weapon physics, early-time fireball growth, and later-time airb1ast phenom- 


ena are calculated with first principle theoretical procedures that combine 


radiation transport phenomena with the hydrodynamic motions of the bomb debris, 


air, and ground media. 
Hydrodynamic calculations have been continued to late 


times and have reproduced fireball rise and toroidal wind flow observed in 


nuclear tests (Ref. 5). 
In principle, such calculational procedures can treat 


the rising dust cloud, the main difficulty being realistic treatment of entrain- 


ment and mixing of air with the fireball or nuclear cloud. 
Credible models for 


ejecta from the crater and of ,dust pickup mechanisms beyond the crater edge are 


also lacking. 
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. Typically, cratering and ground motion calculations are begun by depositing 


some small fraction of the weapon's energy (partitioned between kinetic.and 


internal energy) into a small mass of earth near the point of bur~t. This 


source region is usually described by detailed calculations which consider 


radiation transport as well as shock hydrodynamics. 
The development of a strong 


ground shock begins from this radiation deposition and bomb vapor slap, which is 


largely over within the first microsecond of a surface-burst , nuclear explosion. 


During the first few hundred microseeonds, extremely high pressures exist 


and material resp6nse is modeled as a hydrodynamic fluid (or gas). 
Once the 


shock has progressed into the ground for some distance and the peak shock pres- 


sures have dropped below about 100 kilobars, material strength is included in 


the theoretical model. 
The ground's surface is loaded with an airblast pressure- 


boundary condition~-usually the theoretical airblast loading calculated by 


Brode (Ref. 6) .. Depending on the particular explosion phenomena under study, 


the airblast interaction with the soil may be modeled to account for gas flow 


into cracks, pressure equilibration around rocks and ejecta,airblast-induced 


ground shock-and earth compaction, etc. (Ref. 7). 


It is reasonable t6 question the accuracy of such theoretical procedutes 


since little direct experimental evidence from a high-yield nuclear surface 


burst on rock (or any other surface) is available in these close-in regions. 


There are in fact two basic areas of uncertainty. 
One is the question of numer- 


ical accuracy and the other involves the state of our ignorance of the dynamic 


response of geologic media. 


For example, some studies (Ref. 8) have shown that various first principle 


theoretical procedures treating a given surface burst geometry can produce 


inconsistencies in calculated peak pressures and particle velocities (for 


pressures between 1 and 103 megabars) o~ the order of a factor of 2 or more. 
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For the attenuation rates involved in such diverging ground shock, these incon- 


sistencies imply at least a factor of 2 uncertainty in "equivalent" yield (Ref. 


3) • 


The second major uncertainty concerns the constitutive relations (Ref. 9) 


used to describe the dynamic response of in-situ geologic media, particularly at 


stress levels low enough for material properties to be important, and is compli- 


cated by the difference between the laboratory-determined properties of small 


intact samples and the in-situ response of geologic media. 
Numerous experiments 


in near ..... surface rock media indicate that the dynamic strength of a large rock 


mass is controlled by pre-existing jOints and faults (Ref. 10), and that the 


in-situ-strength may be an order of magnitude less than that· indicated by 


laboratory tests on small intact samples. 


Because of these uncertainties, the most important features of the current 


ground motion calculations are con~idered to be the qualitative features. 
For 


example, such things as the shapes of the peak stress and particle velocity 


contours are probably more accurate than the precise numbers in a given calcula- 


tion. 
Also, a comparison of changes in the results of calculations for slightly 


different physical problems are expected to be more valid than the numbers cal- 


culated in either case. 
In other words, the trends indicated by calculations 


are expected to b.e corre.ct, but quantitative predictions cannot yet .be rigorously 


justified. 


Nevertheless, such improved qualitative understanding can have a significant 


impact on survivability issues even without precise absolute results. 
For 


example, current calculations predict peak stress contours in hard rock that 


differ significantly from early theoretical predictions for a soft rock medium, 


as ShO\V11 in figure 1. 
The more recent results that show higher stress levels 


at a given depth (or conversely, the same stress level. at greater depths) .would 
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Figure 1. 
Comparison of Peak Stress (0.5 kilobar) for 
a Megaton Burst on Rock 


suggest ~hat we should reevaluate our ideas of direct-induced stress wave 


atten~ati6n characteristics in h~rd rock. 
Even this fairly obvious hard rock 


feature, however, cannot be guaranteed on the basis of the calculations alone, 


since no calculations of eArth shock propagation in rock yet include a demon- 


strably realistic model of crack and joint dynamics. 
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SECTION III 


AIRBLAST AND ASSOCIATED PHENOMENA 


In the explosion of a modern l-MT bomb, energy equivalent to 10 15 calories 


is released in a mass on the order of a ton in a time on the order of 10- 8 


seconds. 
This energy is primarily in the form of X rays that radiate from the 


bomb into the air immediately surrounding the weapon, raising its temperature 


to tens of millions of degrees Kelvin and forming the early fireball of hot 


gases. 
The hydrodynamic shock initially created by the bomb vapors overtakes 


the radiation diffusion front at about a millisecond and propagates outward as 


a strong shock wave. 
As the shock moves outward, it diminishes in peak pressure 


and velocity of propagation. 
The uniform pressure in the interior (which at 


high overpressures is 1/2 to 1/3 of the peak pressure)" continues to drop until, 


at about 1.5 seconds, the region interior to the blast wave has dropped to 


ambient pressure and a negative phase begins during which pressures as low as 


3 psi below ambient last for several seconds. 


The peak overpressure (figure 2) passing over or encompassing a hardened 


structure may be largely responsible for the peak accel~rations introduced in 


the structure and, for sufficiently high-frequency structures, may pose the main 


threat. 
However, the motions and deformations imparted to some structures are 


more closely related to overpressure impulse which is defined as the time integral 


of the pressure taken over the duration of the pulse. 
Details of the transient 


airblast environment are discussed in reference 2 and will not be elaborated upon 


here. 


Beyond the range where the hydrodynamic shock temperature drops below SOOOoK, 


air is no longer strongly luminous and the shock front ceases to be an opaque 


barrier to the hotter gases interior to the fireball. At this point the 
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hydrodynamic shock appears to break away from the fireball and the higher 


temperatures of the interior begin to shine through. 
Strong visible light and 


infrared thermal radiation shines on exposed targets for several seconds 


delivering intense thermal energy which results in surface ablation and melting 


of unprotected elements. 
This thermal radiation will continue to shine on the 


target until the fireball gases cool to temperatures too low for effective 


radiant emission. 
If, after the passage of the airblast shock, the target is 


engulfed by the fireball itself, it will be subjected directly to the high 


temperatures that increase behind the shock front as indicated in figure 3. 


These hot gases are by no means stationary and the target may be subjected to 


the strong melting and ablation actions of the hot winds of the.fireball inte- 


rior. 
Boundary-layer effects and local structural geometries may provide some 


protection against these super-hot winds. 


As 'suggested in figure 4, the high winds and drag pressures which follow the 


pa'ssage of the airblast wave present a severe threat to aboveground or protrud- 


ing structures. 
The drag forces may break or deform exposed structures and the 


high winds and entrained soil and rock particles may ablate or destroy sensitive 


surfaces. 
Table I provides estimates of the dynamic pressures, wind velocities, 


temperatures, etc. for several ranges from yields of 0.1, 1, and 10 megatons. 


Of all nuclear environments, the airblast is probably the best understood. 


However, height-of-burst effects are still somewhat uncertain (Refs. 5 and 11). 


Even so, it is believed that near-surface ground motions occurring immediately 


after the passage of theairblast can be predicted with reasonable confidence. 


This level of confidence decreases with increasing time following the initial 


arrival of the airblast-induced ground motions at a given range. 


1. 
AIRBLAST-INDUCED GROUND MOTIONS 


Figure 5 shows the wave fronts in the air-induced ground shock. 
In the 


examples used, the earth media have seismic velocities of 2500 and 5000 ft/sec. 
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and fronts are shown at times when the peak air overpressures are 10,000, 1,000, 


300, and 100 psi. The curves represent positions achieved at uniform seismic 


.. 


velocities and take no account of faster ground-shoc~ propagation at the higher 


stress levels or of variations in seismic velocity with depth. 
It is interesting 


I 


to note that the effect of the rapid slowing of the air shock results in ~ 


steepening of , the wave front at around 300 psi (where it approaches the 5000 


ft/sec seismic speed) with a possible piling up of the signal or waves ~rom a 


c6nsiderab1e range of earlier shock .positions. 
A, similar condition exists at 


about 100 psi for the slower seismic speed example (2500 ft/sec), but is probably 


less pronounced because the air-shock speed decreases more gradually. 
For a 


5000 ft/sec seismic velocity soil and an air shock at 100 psi (traveling at less 


than 3000· ft/sec) , signals in the soil propagate ahead of the air shock producing 


an "outrunning" condition. 


It is obvious, but worth further emphasis, that the wave fronts of figure 5 


do not represent surfaces of ~qua1 pressure. 
In fact, the ,lack of spherical 


"divergence in the wave front directly below the point of burst suggests that 


less geometric attenuation will occur there than in a more symmetric explosion 


(such as a contained burst). 
In the same vein, the ground shock just below the 


shock front at the 300-psi point for the 5000 ft/sec seismic speed cases 


includes signals from pressures considerably higher than 300 psi and could, in 


that region, show ground motions considerably greater than that predicted for 


a simple superseismic air overpressure condition. 


The historically accepted procedures for providing airblast-induced ground 


shock design criteria for protective structures were based on very simple 


intuitive concepts and correlation "of limited data from nuclear experiments 


(Ref 3). 
These procedures have little theoretical basis for vertical moti~ns 


and no theoretical basis for horizontal motions. 'Thus, the application of such 


-empirically based prediction procedures_ must be viewed with skepticism, 
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especially for geologic media other than the soil sites for which t,est data 


exist. 
Because most of the data are for the Nevada Tes~ Site deep dry alluvium 


sites, and because many protective structures are located in quite different 


geologic media, there is a heed for credible prediction procedures that are not 


restricted to the empiricism associated with an inadequate data base. 


In recent years finite difference methods have been used to study airb1ast- 


induced ground motions in layered geologies. 
As suggested in section II, the 


accuracy of these calculations for airb1ast-induced ground motions, is somewhat 


debatable because of questions of numerical accuJ:'acy and also inadequate know1- 


edge of in",:,situ soil and rock behavior under intense stress wave loading. 
How- 


ever, various pa,rametric studies, involving variations in "realistic" material 


properties, have added to our understanding of airb1ast-induced ground motions. 


For example, one parametric study (Ref. 12) that systematica'lly varied the 


geometry and'materia1 properties of a two-layered· elastic-plastic earth model 


indicated that horizontal motions are much more' sensitive than vertical motions 


to variations in yield strength. 
On the other hand, horizontal motions were 


less affected than the vertical motions by the geometric variations ~onsid~red. 


Perhaps the most important result of such studies is the indication that air- 


blast-induced horizontal particle displacements can be considerably larger than 


one-third of the vertical displacements depending on the geology. 
This value, 


which has' often been used as a rule of thumb, may be: satisfactory and even con- 


l 


servative for dry granular soil., However, the re,cent theoretical work suggests 


that unity is a ,more appropriate factor for the more· cormnon soil types with 


water tables near the earth's surface. 


2. 
DEBRIS AND FALLOUT 


Having survived the airblast environment associated with a near miss in a. 


protected installation, one might ask the question: 
What else can happen and 


how soon will, it be over? 
As indicated earlier, immediately following the blast 
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wave positive phase, a negative phase begins in which the winds reverse to blow 


toward ground zero and the overpressure becomes an underpressure (less than 


ambient). 
This negative overpressure can approach as much as 3 psi of suction, 


which could ,exert considerable lift on a sealed, pressurized installation. 


(A 3-psi partial vacuum could list a concrete lid about 3 feet, thick.)' The 


reversed winds may be strong enough to bring back some debris to ciog openings 


or revetments. 


The late fireball is still hot ~ut at nearly normal pressure so that its I 


low-density interior forms a kind of buoyant balloon in the atmosphere. 
This 


several-thousand.;.foot diameter, low-density sphere begins: immediately to rise as 


a bubble as the denser surrounding.air forces it upward. 
The rate of rise 


after a few seconds approaches_400 ft/sec. 
The circulation is such that the air 


velocities in the dust-laden stem that flows up through the rising cloud (figure 


6) are about twice the Cloud-rise velocities, or aSlnuchas 800 ft/sec. 
The 


consequences of such wind velocities can be.better appreciated to,rhen it is'con- 


sidered that the d~ag created by this flow could hold aloft a boulder weighing 


as much as 7 tons or could 10ft lesser rocks and debristo'very high .altitudes. 


The cloud continue~ to ris~ for 4 to 6 minutes, which can take it to altitudes 


'over 60,000 feet, depending on meteorological conditions. 
Even after the cloud 


has stabilized, the stem continues to rise as the circulation persists. 
During 


the time of the initial cloud\rise, much of the cratered debris is aloft on 


various trajectories. 
Much of Fhis material will be excavated at pressures 


below that needed to pulverize or vaporize the rock or soil, and some of H will 


be lofted in essentially its original size and shape. 
If the soil is rocky or 


if concrete and steel structures are involved, some large fragments must be 


expected ~t ranges at least as large as the stem radius; and there is some 


chance that rocks may rain down over a wider area, impacting at near terminal 


velocities. 
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Visibility will be restricted and unpredictable over an area corresponding 


to at least the la-psi distance from such bursts, so that visual assessment of 


the post-burst external environment will not always be possible. 
Direct human 


exposure would be undesirable, possibly even fatal, in the local fallout, which 


(outside the immediate crater area but within 10 miles or so) can rise to 


thousands of roentgens per hour in the first hour, falling to a few hundred at 


the end of a day. 
Total doses (integrated over time) after 18 hours may be in 


exces~ of 3000 roentgens over 1000 square miles. 
Clearly, ~urviving nearby 


surface installations or support structures will not be habitable for many hours 


after a megaton weapon surface burst, even in extreme emergencies. 
The extent 


and intensity of fallout depend critically on weapon design, details of burst 


position, and properties of- soil and surrounding material,as weII' as on the· 


number of bursts in the upwind area . 
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SECTION IV 


FREE-FIELD GROUND MOTIONS FROM SURFACE AND 


NEAR-SURFACE BURSTS 
I 


The ground motions associated with near-surface nuclear bursts are produced 


by the energy coupled directly into the ground at the source point (direct- 


indu_ced ground .motion) and the work done on the ground by. the airblast wave as 


it moves outward over the ground's surface (airblast-induced grouridmotion). 
In 


the past the range to most near-surface structures of interest has been so 


large that only the airblast-induced effects discussed in the previous section 


were considered in atriving at free-field ground motion predictions. 
As near- 


surface structures are designed for more close-in effects, both air-induced and 


, 
direct-induced ground motions must be considered. 
Low-frequency, late-time 


displacements associated with crater formation must also be considered in the 


design ~f shock isblation systems, even at ranges several times the crater 


radius. 


For bursts near the surface of the ground, the airblast loading on the 


ground is essentially independent of the height~of-burst, whereas the direct- 


induced input is quite sensitive to the heightof-~urst. Hence, at a given 


ovcrpreSSUH' 'rnnr,c (l f 'interest, the direct-induced ground motions may or may not 


be larger thl1n [Iir-induced ground motions, depending on the coupling (height ,of 


burst) conditions. 


1. 
INITIAL ENERGY COUPLING 


Only a small fraction of the total yield is effective in producing the 


crater and resulting g~ound motions. 
For example, the kinetic energy requi~ed 


to lift one million tons of ejecta at ioo meters per second (expected from a 


l-MT surface burst) is only about 0 .• 1 percent of the total explosive energy of 
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the weapon. 
Ground motions at late times also involve quite small percentages 


of the explosive energy. 
Thus, the details of how energy is-coupled into the 


ground can be very important to the environment relevant to protective structures 


even though they affect'the overall energetics in a minor way. 


~ 
If the bomb explodes within a meter or so 'of the surface of the earth, then 


the downward-directed portion of the X':ray yield contributes a short pulse of 


high-energy photons to the ground surface just before the, expanding bomb debris 


arrives. 
In a fraction of a microsecond, the reem:1ssion of radiation from the 


surface of the ground and effective penetration of the radiative diffusion front 


into the: ground are essentially comp+eted. 
Most of this radiatively coupled 


energy will be distributed in a thin tapered disc of material a few centimeters 


thick. 
The exact depth depends on the radiative characteristics of the bomb and 


height of- burst. 


The X-ray-heated ground material vaporizes and blows up. 
Some of it collides 


with the expanding weapon vapors. 
The interaction of the bomb debris with earth 


blow-off .tends to contain the blow-off and to sustain the shock induced in the 


ground by the X-ray heating. 
After about a microsecond the radiative diffusion 


into the ground is no longer very important and the effective radiative energy 


for producing ground shock has been brgely converted into internal energy in 


the shallow source region near'thesurface of the ground. 
Thus, a source of 


internal energy (due primarily to radiative coupling) and kinetic energy 


(principally from bomb dehris) may be defined as i~put for subsequent hydrody- 


namic calculations. of the ground shock. 
Figure 7 shows typical initial condi- 


tionsfor the ground motion problem. 


, The initial partitioning of energy between kinetic and in~ernal energy is 


important because downward-directed kinetic energy from the bomb debris is more 


effective in producing ground shock than the isotropic deposition of the same 
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1 MT 
-- 1 p,sec 


HOl air: 


Bomb: 


Radius -- 34 m 


Temperature - 
2.000.000oK 


Pressure - 
130 kb 


Velocity - 
I kmlsec 


t 


Radius - 
2 m 
Temperature - 
4.000.000oK 
Pres5ure- 
lo~'kb 


Velocity - 
1000 kill/sec 


Figure 7. 
Initial Conditions for Surface-Burst Cratering Motion 


amount of internal energy. 
For example, a downward-directed kinetic energy 


source produces close-in pressures that are about twice as high as those pro- 


duced by the same energy all in internal energy in the same mass source (Ref. 8). 


The height of the burst has an important effect on the radiative coupling. 


The air acts as an energy sink for out'ward-directed photons, converting them to 


thermal energy in a small firebali about the burst point. 
Most of this early 


fireball energy .is reemitt~d in relatively low-energy photons, but on a much 


longer time scale than that usually associated with the earth coupling phenomena. 


For heights of bursts greater than about 100 feet the prompt ~adiative coupling 


. , 


for most weapons is expec.ted to be negligible. 
For a 50-foot height of burst 


the initial radiative coupling is expected to be on the order of one-half that 


for a contact burst, and the kinetic energy in the debris will be largely 


dissipated before reaching thejground. 
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2. 
DIRECT-INDUCED GROUND MOTIONS 


As an indication of the importance of shallow_burial on hydrodynamic coup~ing, 


consider the results from hydrodynamic calculations for a 790-KT source (without 


radiation coupling considered) at 0, 1 meter, and 8 meters depths of burial 


(Ref. 13). 
As is seen in figures 8 through 10, the depths and ranges near the- 


surface reached by given pressure levels increase with depth of burial. 
Note 


the drastic increase in the "pressure" near the free surface due to just an 


8-meter depth of burial--about a factor of 7 at 40 meters range. 


In addition to the above-mentioned hydrodynamic coupling effects, slight 


penetration of a nuclear weapon into the ground's surface (perhaps as little as 


3 to 5 meters for a megaton burst) has an even more dramatic effect on the 


radiative coupling. 
The radiative coupling may be increased by more 'than a 


factor of 5 in going from a contact burst to one with shallow burial. 


To estimate deep undergrounds motions, designers and engineers, .have used the 


concept of an effective or equivalent yield (W ) defined as being the yield of 
- 
-, 
. 
e 


a contained nuclear explosion that would provide the observed peak range 


beneath ground zero of a nuclear surface burst of yield W. 


then defined as an effective coupling factor. 


The ratio W /W is 
e 


, This useful design concept should ,not be confused with the actual energy 


coupled into the ground at early times. 
Although it is expected that ~here may 


be some relationship between the two, that relationship is neither constant nor 


well understood. 
All that can be said at present is that if the initially 


coupled yield is increased, the equivalent yield is also increased. 
In any 


event, the equivalent yield should not be construed to imply anything about the 


energy delivered to the ground. 


Peak-particle velocity measurements fro'm nuclear explosions buried in hard 


rock (granite) and soft rock (tuff) are shown in figure 11 (Ref. l4)~ A best 
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Peak Pressure Isobars for Surface 790-KT Source 
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Figure 9. 
Peak Pressure Isobars for I-Meter Burial of a 790-KT Source 
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Figure 10.- 
Peak Pressure Isobars for 8-Meter Burial of a 790-KT Source 
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. estimate particle velocity (v) 
i~ 


[ 
32 


R 


8, We 1/3J2 
v ~ 1.lC 
1000 


fps 


where R is the range (ft)' and C is the seismic compressional wave speed (fps). 
, 


Based on peak-particle velocity data from low-yield nuclear detonations in 


underground cavities in granite and tuff, an equivalent yield coupling factor of 


about 10 percent would be inferred. 
This measured equivalent yield coupling 


, 
factor is about twice the recommended value for large yield nuclear devices 


,'>. 
'\ 
• 


(Ref. 3). 
The more massive (low yield to mass ratio) low-yield devices are 


expected to more,efficiently couple,energy -to the ground than the higher energy_ 


density, large yield devices. 
Thus, we will assume an equivalent coupling 


factor of 5 percent consistent with previous estimates, and estimate the ground 


motions beneath a surface burst on a rock medium. 


'The free surface perturbs the motion field to, the extent that a spherical 


contour assumption becomes progressively worse· as the free surface is approached. 


As noted in figures 8 through 10, hydrodynamic calculations suggest·that the 


strong shock produces an early time spherical stress and motion field out to 


about 60 degrees from the axis of symmetry. 
Beyond about 60 degrees,the free- 


surface relief destroys the spherical symmetry and leads to ,!owerpeak stresses 


above that angle. 
For pressures less than/about 100 kilobars, soil and rock 


strength propert.lcH are expected to be increasingly important in modifying even 


the near axi~. ~:ph'~ric9l field. .Although the shock front appears to be spherical 


(in a homogeneous medium), the combination of free-surface and strength effects 


lead to a reduction in peak stress as the field point on the, shock front moves 


away from the axis,of symmetry. 
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At lower stress levels, the field of motion below a surface burst might be 


thought of a consisting of three regions (Ref. 15), as shown in figure 12. 


Region I consists of a spherical field and extends to a conical surface with an 


apex half-angle of about 30 degrees. 
In region II (between about 30 and 65 


degrees), tqe early-time principal stress directions are as in a spherical 


field, but the !.'!tress amplitudes are significantly different than in a spherical 


field. 
Above about 65 degrees the air-induced surface motion is coupled with 


the direct-induced motion so that even qualitative statements- are difficult. 


v 
-DffiECT 
I INDUCED MOTION 
I 


REGION II 


REGION I 


·Affi INDUCED 


MOTION 


REGION m 


Figure 12. 
Phenomenological Regions for the Surface-Burst Geometry 
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Based on several theoretical studies for hard rock and some simple concepts 


from geometrical optics theory, it is assumed that 'direct-induced peak-particle 


velocity (v) and stress (a) contours can be represented by 


v~ l.le 


1000 


(J = pcv 


(cos e + 0.1 sin 8) (fps) 


where Rand e are the usual spherical coordinates and W is the equivalent yield 
e 
in KT. 
This fit is consistent with theoretical calculatipns!of surface bursts 


on rock media. 
The direct-induced peak velocity contour shape given by 


F(e) ='cos e + 0.1 sin e has a slight bulge at 8 = 5.7 
0 where F(e) = l.09. 
Thus, 


Fee) does not have precisely the correct behavior and might be thought of as the 


leading" term approximat:j.on in a series expansion of the "correct". function 


representing the contour. 
It should be noted that the airblast-induced ground 


motion alters the peak values near the earth's surface, leading to contours with 


the shape sho~~ in figure 1. 


'The prediction of peak-particle displacement and other late-time phenomena 


is a considerably more, formidable task than predicting peak-particle velocities, 


stresses) and strains. 
Generally, these latt~r variables, are controlled by the 


inunediate response of the earth media to the shock-wave passage; whereas, peak 


displacements are governed by late-time behavior of the in-situ rock or soil. 


This late-time phenomena can be affected by reflections from even rather distant 


,inhomogeneities and by the in-situ behavior of the earth media. 
The close-in 


late-time motions in rock are strongly affected by the motion along jointing 


surfaces. 
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Note that peak-particle velocity data from underground experiments in Ranier 


Mesa tuff (a "soft" rock) and granite (a "hard" rock) (shown in figure 11) differ 


by a factor of about 3.5 (roughly consistent with the ratio of seismic velocities). 


This difference is consistent with the order of magnitude difference in confined 


modulus (pc 2) one would expect based on laboratory tests of small samples of 


granite and tuff. 
However, as shown in figure 13, an order of magnitude differ- 


ence between the displacements in tuff and granite is not indicated--a fact that 


1 


would be surprising if one expected in-situ rock behavior to be that based on 


laboratory tests of sm9ll samples'. 


It should also be noted that the cavity displacement in French underground 


tests in granite were significantly smaller than those in experiments in Nevada 


Test Site granite (Ref. 16). 
Some have attributed this difference to the differ- 


ence in water content--the French tests being in much dryer rock. ' In any case, 


r 
' 
there is the suggestion that displacements in granite in par~ of thiw6rld are 


different from displacements in granite at other parts of the world .(and further, 


that displacements in. some granite and some tuff are comparable). 


The previOus discussion involved the manipulation of peak-particle velocity 


data from contained bursts with little regard to free-surface effects (excep~ 


for peak stress and particle-velocity contours such as those in figures 8 through 


10). 
To demonstrate late-time mOtion effects, we will consider some typical 


data from ·HINERAL ROCK,a lOa-ton high-explosive spherical charge that was 0.1 


of the charge radius buried (Ref. 17) in granite. 
This experimental source 


geometry was chosen in an attempt to partition the airblast-induced and direct- 


induced ground motion in a proportion like that expected for a nuclear surface- 


burst. 


, 


At a horizontal range of about 2 crater.radii (==70 feet), the peak-particle 


velocity data scatters between about 3 and 15 fps. 
Typical particle velocity 


time histories at a 70-foot range are shown in figure 14. 
The upward and outward 
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Figure 13. 
Comparison of Peak-Particle Displacement from 
Contained Bursts in Hard and Soft Rock 
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Figure 14. 
Mineral Rock Particle Velocity at a 70-Foot Range 
and a 2-Foot Depth 


direct-induced ground motions completely dominate the ground motion. 
Note that 


there is little "elastic" recovery indicated by these data traces. 
While 


q~estions of base line shifts always cause problems in the analysis 6f integrated 


transient accelerometer and velocity·gagedata, there can be little doubt of 


predominantly upward and outward particle-displacements. 
The majority of per- 


manent displacement data obtained from pre- and post~test surveys are shown in 


figure 15. 


Based qn such information as that discussed above, the near-surface transient 


direct-induced ground motions are predominantly upward and outward. 
Peak ver- 


tical (upward) particle velocities are expected to be somewhat less than peak 


horizontal particle velocities, but peak vertical (upward) and horizontal par- 


ticle displacements are expected to be about equal. 
Considerable data scatter 
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Figure 15. 
Free...,Surface Permanent Displacements (Mineral Rock)· 
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about nominal or best estimate values should be expected (factor of 3 in par- 


ticle velocity and factor of 4 in particle displacements). 


3. 
CRATER-INDUCED GROUND MOTIONS 


Ground motions associated with the relatively late-time crater forming pro- 


cesses are clearly large at the edge of the crater and can conceivably dominate 


the late-time surface motions at distances large with respect to the crater 


radius. 
Evidence that such might be the case was the observation of frozen 


surface waves in a radial asphalt strip at about 2 crater radii from Distant 


Plain 6 (Ref. 18), a 100-ton HE experiment at the Suffield Experimental Station 


(figure 16). 


Based o~ the notion th~t such motions should be relatable to some character- 


istics of the crater, we have analyzed the maximum transient and permanent 


displacement data from a large number of high-explosive surface bursts. 
Based 


on this analysis, we believe that transient displacement data from all above- 


surface crate ring bursts can be compressed into a single scaling curve by use of 


a length scaling factor equal to the cube-root .of the apparent crater volume. 


A best estimate fit to this data is 


= 0.45 V4 / 3 (ft) 
R3 


where ~ and dv are peak transient horizontal (outward) and vertical (upward) 


displacements, V is the apparent crater volume, and R is the horizontal range 


from ground zero. 
Assuming V ~3x 10 7 ft 3 per megaton for hard rock, the above 


formula predicts displacements as shown in figure 17. 


The displacement pulse under discussion has quite a low-frequency content. 


For example, note·' the,~ypical pulses shown in figure 14. 
For such 100-ton high- 


explosive experiments in rock, the rise time of the displacement pulse is on 
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Figure l7~ Estimates of Crater-Induced Displacements 
from Contact Bursts on Hard Rock 


the order of 0.1 second and in soil the rise time is the better part of a 


second. 
For nuclear bursts it is expected that the horizontal displacement 


rise times might be even longer--on the'order of a few 'seconds. 
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SECTION-V 


SUMMARY AND CONCLUDING COMMENTS 


This review of blast and shock related close-in nuclear weapons effects has 


emphasized recently obtained information not generally available, and it comp1e- 


ments the general surveys given in references 1 and 2. 
Perhaps the most impor-' 


tant recent development is the application of first-principie theore)tica1 pro- 


cedures to study close-in ~uc1ear phenomenology. 
In particular, such theoretical 


results are increasing our understanding of airb1ast and ground-shock phenomena; 


thereby providing a bridge for wider application of the limited nuclear test 


data base. 


These studies lead us to question some of the early notions that resulted 


from an incomplete understanding of basic phenomena and the incorrect extrapo1a- 


tion of Nevada Test Site ground-motion test data. 
For example, horizontal air- 


blast-induced ground motions in some geologies are likely to be larger than 


est·imated in the early 1960s.·· The continuing interest in designing survivable 


_ systems close-in to surface bursts has led to, theoretical and experimental 


programs to improve our understanding of cratering and associated,c1ose-in 


phenomena. 
Although the first comprehensive cratering calculation was performed 


over 10 years ago, c~rrent theoretical studies still provide only qualitative 


answers and must be ~omp1ementea with well-conceived experiments. 


In partlculLlr, extl'llt;( ve experiments have been conJucted to simulate the 


I 
.close-in ground-shock phenomena in hard rock. 
The key lessons learned from this 


program are that in-situ near-surface hard rock responds to intense stress-wave 


loading in a manner quite different than would be expected on the basis of many 


previous theoretical studies of highly idealized homogeneous materials. 
The 


real world contains fractures and joints that dominate the dynamic response 
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characteristics of large rock masses. 
Relative motions occurring along such 


pre-existing planes of weakness produce spatial discontinuities in the displace- 


,ment field that can pose a severe problem for hardened structures, particularly 


if the designer did not understand the environment. 
Although this phenomenon 


is not predictable in a precise sense, a substantial data base provides high 


confidence that it will be produced by intense stress waves in many locations. 


(Much of the data scatter indicated in figures ll, 13, and 15 is probably caused 


by such a nonuniform spatial motion field.) 
Thus, ground motions produced by 


high-yield nuclear contact bursts on hard rock are expected to vary azimuthally 


in a manner that might possibly be associated with the pre-existing joint pattern, 


assuming ground zero is determined. 
Unfortunately for the analyst the location 


of a threat burst point is inherently uncertain; thus such predictions are only 


of academic interest. 


Thus, the most significant conclusion is that a·considerabJ.e uilcertainty 


will always exist in providing the' criteria for the design of hardened stru.ctures. 


Some uncertainty is associated with our inability to model what is known and may 


be reduced by additional research efforts. 
However, an inherent uncertainty is 


associated with our inability to ever know all important geologic details (such 


as local inhomogeneities) that produce the observed data scatter in figures 11, 


13, and 15. 
The design of facilities to survive such close-in environments must 


recognize this uncertainty and attempt to use procedures that circumvent their 


possible destructive effects. 
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